Bragg solitons in nonlinear PT-symmetric periodic potentials 
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It is shown that slow Bragg soliton solutions are possible in nonlinear complex parity-time (VT) 
symmetric periodic structures. Analysis indicates that the PT-symmetric component of the periodic 
optical refractive index can modify the grating band structure and hence the effective coupling 
between the forward and backward waves. Starting from a classical modified massive Thirring 
model, solitary wave solutions are obtained in closed form. The basic properties of these slow 
solitary waves and their dependence on their respective 'PT-symmetric gain/loss profile are then 
explored via numerical simulations. 
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I. INTRODUCTION 

Periodic structures play an important role in the 
general area of optics. As in solid state physics, the peri- 
odicity in their refractive index can lead to a succession 
of photonic band gaps and transmission bands. In many 
applications, these properties are used to obtain high 
reflectivities, frequency filtering, and high-dispersion 
characteristics [l|. Index gratings, whether in the bulk 
or embedded in optical fibers d, 0| , are examples of 
such structures. Even though in principle these periodic 
configurations can always be rigorously analyzed using a 
Floquet-Bloch approach, on many occasions a coupled- 
mode formalism will suffice. As shown by Kogelnik, 
this latter formalism is particularly successful when 
the periodic index perturbation is weak, in which case 
the coupling between the forward and backward waves 
occurs over a distance of several wavelengths Q. In this 
regime, the interaction can be described through the 
so-called slowly varying approximation, which in turn 
leads to a relatively simple system of coupled equations. 

The behavior of Kerr nonlinear optical periodic sys- 
tems was first addressed in 1979 under continuous-wave 
conditions in conjunction with optical bistability @. 
Few years later, it was realized that this same system 
can also support a special class of soliton solutions-the 
so called Bragg solitons 041l- Unlike optical solitons 
propagating in nonlinear dispersive fibers, this family 
of waves is made possible by nonlinearly interlocking 
both the forward and backward propagating modes 
fiol [Tl| . In doing so, these wavepackets "open up" a 
defect band within the forbidden band gap thus allowing 
energy transport. Given that under Bragg conditions 
this propagation is linearly forbidden (the grating acts 
like a distributed mirror), the resulting propagation can 
be very slow as in the case of self-induced transparency. 
In general, the velocity of Bragg solitons can range 
from zero (fully immobile light) to c/n, depending on 



excitation conditions. We would like to emphasize that 
so far, this class of waves has been primarily inves- 
tigated in conservative systems. The question arises, 
as to how they will behave in non-conservative envi- 
ronments, especially in the presence of linear gain or loss. 

Parity-time (VT) symmetry in op tics has recently 
attracted considerable attention [12 29]. While VT 
symmetry was first explored within the quantum domain 
|30l - [33l ] , it is in optics that has found a straightforward re- 
alization where its implications can be directly observed 

and studied 00. As sh ° wn 

in 0, [l3j, an optical 
system obeys VT symmetry provided that its complex 
refractive index distribution n(r) = n^(r) + m/(r) 
satisfies the condition n*(r) = n(—r) . In other words, 
the real index profile must be an even function of 
position while the gain/loss must be odd. It can be 
shown that for such structures, a real propagation 
constant (eigenenergies in the Hamiltonian language) 
exists for some range (exact VT- symmetric phase) of the 
gain/loss coefficient. For larger values of this coefficient 
the system undergoes a spontaneous symmetry breaking, 
corresponding to a transition from real to complex spec- 
tra (broken VT -symmetric phase). The phase transition 
point, shows all the characteristics of an exceptional 
point (EP) singularity. Abrupt PT-symmetry breaking 
has been recently observed in both active and passive 
experimental arrangements fl4l . [l5l |28| . In addition non- 
reciprocity in propagation as well as double refraction 
and energy oscillations have been predicted in periodic 
lattices and coupled structures. The possibility of 
unidirectional invisibility was put forward in linear and 
nonlinear PT-symmetric gratings and the properties of 
PT-symmetric scatterers and lasers were also discussed 
[23l 0). Finally, along somewhat different lines, the 
prospect for static optical solitons in PT-periodic arrays 
has been considered in several studies |34|-|38[ . 

In this work we demonstrate that a new family of 
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optical Bragg solitons is possible in Kerr nonlinear 
PT-symmetric periodic structures. Star ting from a 
classical modified massive Thirring model [331, solitary 
wave solutions are obtained in closed form. The basic 
properties of these slow solitary waves and their de- 
pendence on their respective PT-symmetric gain/loss 
profile are explored and pertinent numerical simulations 
are carried out to elucidate their behavior. We also show 
that at the exceptional point, the evolution equations 
decouple, thus allowing a special class of solutions. 



II. THEORETICAL ANALYSIS 

We begin our work by considering a PT-symmetric 
optical grating having a periodic complex refractive index 
distribution: let us consider a fiber with the following 
refractive index of the core: 

n = no + niRCOs( — z) + inusin{—z) + n 2 \E\ (1) 

In this profile the first term stands for the refractive index 
of the background material involved while the three other 
terms are considered to be small perturbations on no; the 
second term describes periodic Bragg grating, the third 
term represents the superimposed complex VT potential 
(gain or loss) and the last term accounts for the Kerr 
nonlincarity. We now express the solution as a sum of 
forward and backward propagating waves: 

E = Ef(z,t)exp[i(f3oz — ^ot)]+E b (z,t)exp[—i(f3 z+ujot)] 

(2) 

where luq = 2ttc/\q is the carrier angular frequency, Ao is 
the free space wavelength and (3o = ^o w o/ c is the un- 
perturbed propagation constant. Finally Ef(z,t) and 
E b (z 1 t) represent slowly varying amplitudes for the for- 
ward and backward waves respectively. In this case, it 
can be directly shown that the two slowly varying enve- 
lope functions satisfy the following coupled wave equa- 
tions: 



JdEf IdEA _ i2S 
\ dz v dt J (3a) 
+ 1 (\E f \ 2 + 2\E b \ 2 )E f =0, 

V dz v dt J V yj 5 (3b) 

+ 7 (|^| 2 + 2|£; / | 2 )£; b = o, 

In the above equations v = c/no is the wave velocity in 
the background material, k = tttiir/Xo is the coupling 
coefficient arising from the real Bragg grating itself, and 
g = irnu/Xo is the anti-symmetric coupling coefficient 
arising from complex VT potential term. In addition, 
8 = [tiq/c){u)q — u>b) is a measure of detuning from the 
Bragg angular frequency ujb = 7rc/(noA) and 7 = n2Wo/c 



is the self-phase modulation constant. 
In the linear regime, the properties of Eq. ([3]) can be 
readily understood by using the following gauge trans- 
formation, Ef = Fe~ lSz e iv5t , E b = Be lSz e ivSt , in which 
case one obtains: 

+ *(!! + iw) + (* + .9)5 = 0, (4a) 
-i(W-m+(^-9)F = 0, (4b) 

By assuming time harmonic solutions of the form, 
{F,B) = (Fo, Bo)exp (i(Kz — fit)) we arrive at the dis- 
persion relation: 

A' 2 = ^-(^ 2 -.9 2 ) (5) 
v 

The effect of the PT-symmetric term arising from g on 
the overall dispersion characteristics of this Bragg grat- 
ing is obvious. In essence, its presence can effectively 
shift the photonic band gap as illustrated in Fig. [TJ for 
different ratios of g/ 'n. In Fig. [T] the dispersion prop- 




FIG. 1. (Color online) Band structure of a PT-symmetric 
periodic grating (linear case) for different ratios of g/n; (a) 0, 
(b) 0.8, (c) 1, (d) 1.2 



erties of this periodic PT grating are depicted for three 
different regimes, depending on the ratio of gjn\ b) for 
g < k (below PT-symmctry breaking threshold) the 
band structure has essentially the shape of an ordinary 
Bragg grating-with the photonic band gap reduced, c) 
for g = k (at the VT threshold or exceptional point) 
the band gap is closed and the dispersion curve is identi- 
cal to that expected from the homogeneous background 
material, and d) for g > n (above threshold) where no 
band gap exists and the dispersion relation is totally dif- 
ferent in shape. As Fig. 1(d) illustrates, above the VT- 
symmetry breaking threshold, around the origin, there 
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is always a range of wavevectors associated with com- 
plex frequencies. As we will see, this latter observa- 
tion explains why in this case field configurations can 
grow/decay exponentially with propagation distance. In 
addition, in this same regime the group velocity is always 
larger than velocity of light within the background ma- 
terial. In this work, we mainly restrict our attention in 
the first range, i.e., we will assume that the VT grating 
will be operated below the VT threshold where the entire 
frequency spectrum is real. 



III. NONLINEAR DYNAMICS AND SOLITARY 
WAVE SOLUTIONS 

In this section we investigate the existence of solitary 
wave solutions for the coupled wave Eqs. [3j To do so, 
we exploit the existing similarity between Eqs. [3] and 
of that of the massive Thirring model [3!|. By intro- 
ducing the two parameters p = y/ (k — g)/{n + g) and 
= \J k 2 — g 2 and by employing the gauge transforma- 

pBe i5z e ivSt , these coupled 



tions E f = Fe~ l5z e lvSt , E b 

wave equations can be written in the following form: 

+ i (t? + ~lr) + + 7 (\F\ 2 + 2p 2 |B| 2 ) F = 0, (6a) 
- i (If - ) + «rF + 7 (P 2 \B\ 2 + 2\Ff) B = 0. (6b) 

We note that the above mentioned gauge transformation 
is only valid when k > g, e.g. below the VT threshold 
point. As a next step we consider a solution of the form: 



(7) 



where the constant a and the function n(z,t) remain to 
be determined. On the other hand, (ij'f^b) represent 
solutions to the Thirring model |8HlO(: 



(8b) 



ipf = + J^±sin(cr)e l ' s, sech (9 - if) 



tpb = — ./ ^ £ -Asin(cr)e l<I> sec/i I 



where <f> and 9 are functions of z and t defined as follows: 

z — vmt 



9 = Kpsin(a) 



$ = K p cos(a) 



v/l-m 2 



vt 



(9) 



(10) 



In the above, the dimensionlcss quantity m is defined as 
m = (1 - A 4 )/(l + A 4 ) and finally A and a (0 < a < tt) 
are free parameters. After inserting these solutions into 
Eq. [H]we then obtain: 



dr) 
drj 



I (la 2 | „2„2 
= + {2 AT +P 



I) sin(a) \sech(0 - if )| 2 ,(lla) 
(|q 2 P 2 A 4 + a 2 - 1) sin(a) \sech{9 - if )| 2 (llb) 



A valid solution of Eqs. [IT] requires that both sides are 
equal. This condition in turn determines the unknown 
coefficient a: 



1 



1 



4A 4 



-1/2 



(12) 



Finally rj can then be obtained by integrating either one 
of Eqs. [TTJ 

77(0) = 2 + p 2 a 2 - lj tan^ 1 (tanh(6)tan 

(13) 

Here it is worth discussing the velocity and instantaneous 
frequency associated with this soliton solution. Accord- 
ing to Eqs. ((HJ and ((!]) the soliton velocity can be readily 
obtained from: 



1 - A 4 
1 + A 4 



(14) 



Hence the soliton velocity can reach any value between 
zero (A = 1) and the group velocity in the background 
medium (A = 0). Using an amplitude and phase repre- 
sentation of Eqs. and flSJ), the corresponding phase of 
this soliton solutions could be written as, 



a = 77 ■ 



$ ± tan 1 (tanh(9)tan ( — 



(15) 



where the plus and minus signs correspond to the for- 
ward F and backward component B respectively. Note 
that these phases are obtained after the aforementioned 
gauge transformation. Hence to obtain the actual phases 
for the forward and backward waves (Ef,Eb) the term 
vSt^fdz must be added to these phases respectively. The 
instantaneous angular frequency can then be obtained 
from a first order term Taylor series expansion of the 
respective phase of Eq. (fl~5j) : 



VI - m 2 

K,pV 

VI - m 2 



cos(a) 



2A 



1 2m tan (%) sech 2 <8<z, t = 0)) 

+ p 2 a 2 -l±0.5 



1 + tan 2 (f ) tanh 2 (8(z, t = 0)) 



(16) 



Given that a gauge transformation was used, the quantity 
vS must be subtracted from the result of Eq. (fT6| , which 
is measured with respect to the carrier frequency. Thus 
the total instantaneous angular frequency of this soliton 
solution is given by lo s = Q s —v5+luq = Q s +lub- Accord- 
ing to the linear dispersion analysis used in the previous 
section, the frequency band gap for the "PT-symmetric 
grating can be obtained from —k p v < Q < k p v. There- 
fore, based on Eq. (|16p the soliton frequency fl s may or 
may not lie in the band gap. 

Up to this point, the solutions were obtained for k > <?, 
i.e. before the VT symmetry is broken. On the other 
hand, at exactly the 'PT-symmetry breaking point [k = 
g), the effective coupling coefficient k p goes to zero. In 
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this case, the evolution equations are not completely de- 
coupled and can be more effectively treated in the orig- 
inal set of variables. By introducing the gauge transfor- 



mations E 
wave Eqs 



= F e _taz e 1 
reduce to: 

1 dF 



uSt 



E h = Be 



i8z „ivdt 



, the coupled 



+ i (H + £§D + 2 * B + 7 (l^l 2 + 2|B| 2 ) F = 0, (17a) 
-*(H-^lf)+7(|B| 2 + 2|F| 2 )B = 0, (17b) 

The linear coupling term between the forward and back- 
ward waves now breaks the symmetry in the evolution 
equations. Note that there is no energy transfer from the 
forward wave to the backward but the backward wave 
facilitates energy transfer to the forward. This can be 
better understood by considering the general solution of 
Eq. I17b[ given by: 



B = b(y)exp ( —17 



b 2 (y)x 



\F\ 2 d£ 



(18) 



where 



Uz - vt), 



y 



h(z + vt) are forward and 



backward propagation coordinates and b is an arbitrary 
function. On the other hand Eqs. 1 171 admit a trivial solu- 
tion when B = 0. In this latter case, Eg. I17al reduces to 
that describing a forward propagating wave in the pres- 
ence of nonlinear self-phase modulation, which admits 
the following solution: 



F = f(x)exp(i-yf 2 (x)y) 



(19) 



where / is an arbitrary function. In the other words, in 
this regime the intensity profile of the forward propagat- 
ing wave remains invariant during propagation while no 
energy is transferred to the backward mode. 



IV. NUMERICAL RESULTS 




FIG. 2. (Color online) Propagation dynamics of a solitary 
wave solution in a PT-symmetric Bragg structure; intensity 
evolution for both the forward (left) and backward waves 
(right) during propagation. 



solitons is taken to be ir/2, and parameter A that deter- 
mines the common velocity of the two constituent waves 
is taken to be 0.8. Figures (|3l4p on the other hand show 



Forward 



Backward 




In this section we exemplify our results through nu- 
merical simulations of Eqs. ([6]). The numerical methods 
used for solving the coupled wave equations presented 
are based on finite difference methods using different 
discretizing approaches in order to account for numerical 
stability [3, [l0|. Here for discretization we use Euler's 
method that is based on a first order approximation 
for both temporal and spatial derivatives. In this case 
stability would not be an issue as long as the temporal 
step size is way smaller than the spatial step size. First 
we investigate the behavior of the solitary wave solution 
given by Eqs. ((7] fTS|) . Figure [2] depicts the corresponding 
propagation dynamics of this solution for both the 
forward and backward waves. According to this figure, 
these two components propagate at a common velocity 
and they have the same profile (except from a scaling 
factor that is clear from Eqs. ©). In this numerical 
example g/n = 0.8, and the space-time coordinates are 
normalized as follows: Z = kz and T = nvt. In addition 
the forward and backward electric fields arc also here 
normalized with respect the quantity Eq = yfizpy. The 
parameter a that determines the beam width of these 



FIG. 3. (Color online) Propagation dynamics of a Gaussian 
wavepacket when injected only in the forward direction when 
the VT grating is operated below the PT-symmetry breaking 
threshold. (a),(b) depict the forward and backward compo- 
nents respectively and (c) the associated energy as a function 
of normalized time. 

the evolution of a Gaussian pulse when it excites only 
the forward wave within such a 'PT-symmetric Bragg 
grating, for two different cases: below the PT-symmetry 
breaking point and at threshold. In these simulations 
g/n is set to be 0.8, 1 respectively. In these figures the 
total energy of each component that is proportional 
to \H(z, t)\ 2 dz (where H is either a forward or a 
backward wave) is also plotted as a function of time. In 
the case of PT-symmetric soliton solutions this quantity 
is constant with propagation. 

According to Fig. [31 below PT threshold there is an 
oscillatory power exchange between the forward and 
backward waves. In this same regime, by increasing 
the amplitude of the imaginary potential (amplitude of 
gain or loss), the rate of this energy exchange decreases. 



5 




5 10 

time 



FIG. 4. (Color online) The same as figure [3] when the VT 
grating is operated at the 'PT-symmetry breaking threshold. 
(a),(b) depict the forward and backward components respec- 
tively and (c) the associated energy as a function of normal- 
ized time. 

Figure |4] on the other hand, shows that the forward 
Gaussian wave remains unchanged during propagation 
while the backward wave is not excited at all. This 
is in agreement with our previous discussion, as ex- 
pected from equation (1191) . This is because there is no 



energy coupling between the forward and backward wave. 



V. CONCLUSIONS 

In this work we have demonstrated that a new family 
of optical Bragg solitons is possible in Kerr nonlincar 
■pT-symmetric periodic structures. By considering the 
connection to the classical modified massive Thirring 
model, solitary wave solutions were obtained in closed 
form. The basic properties of these slow solitary waves 
and their dependence on their respective "PT-symmetric 
gain/loss profile were explored and pertinent numerical 
simulations were carried out to elucidate their behavior. 
Finally, of interest will be to examine if similar concepts 
can be applied in other periodic structures as for 
example in nonlinear optical mesh lattices [13, EJ . 
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